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1. Introduction 

Lattice diagrammatic techniques can be valuable tools to obtain insight into the strong cou¬ 
pling limit of QCD and related theories. We consider a particular diagrammatic approach which 
was introduced in the 80’s to study chiral symmetry breaking in QCD at infinite coupling, as 
Nf 0, in [1], and then further developed in [2], More recently this approach has been picked 
up again to address the question of chiral symmetry restoration in the case of QCD with a large 
number of fermion flavours Nf. In particular, the simulation results in [3] for the chiral condensate 
at infinite coupling as a function of Nf show evidence of a first order transition to a chiral symme¬ 
try restored phase at a critical value of Nf ~ 13 staggered flavours. Although such a transition is 
well documented at more moderate coupling strengths, its presence at infinite coupling came as a 
surprise, because analytical calculations based on a 1 /d expansion [4], or mean field [5], suggested 
that chiral symmetry would remain broken for all Nf at infinite coupling. The lattice diagram¬ 
matic technique of [1, 2] was then reintroduced and extended to account for contributions arising at 
nonzero Nf in [6]. There are two solutions for the normalised chiral condensate as a function of Nf 
obtained in [6]. One of these solutions matches onto [2] in the Nf — > 0 limit, where the normalised 
chiral condensate goes to ~ 0.66 as Nf —* 0, then increases in magnitude as Nf increases. The 
other solution goes to infinity as Nf — > 0, and decreases as a function of Nf. For both solutions, 
there is a common critical value of Nf ~ 10.7, beyond which only complex solutions for the chiral 
condensate exist. It would be good to understand this better. The idea of this note, and of our recent 
longer paper in [7], is to develop a procedure, inspired by [6], which can be used to calculate the 
chiral condensate by collecting the contributions from all possible diagrams which can be formed 
out of a truncated number of sub-diagram types. 


2. (t/At/r) at g — oo 


As in [6], we begin by generalising the procedure in [2] to incorporate contributions which 
arise at nonzero Nf. The iterative procedure we employ to generalise [2] is different from that of 
[6], and we summarise it below using the notation of [2, 6]. 

Integrating out the fermion fields puts the chiral condensate in the form 


(V'MV'to) = - lim tr 

m-> 0 


f dU det [l + K~ l M(U)] [l + K~ 1 M(U)] l K 


JdU det [1 +K~ l M(U)\ 


(2.1) 


with 
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( 2 . 2 ) 


for /i = 1 The form of (i/ai/a) in (2.1) suggests expanding in powers of K 1 M, resulting in 
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The presence of the trace in (2.3) and (2.1) allows for simplifications using tr [odd # of y+s] = 0, 
such that only contributions from terms with ( K~ l M) n for n even are nonzero. In addition, due to 
the SU (N c ) integrals over the U’s, the only nonzero diagrams are those where each link has U m U ' n , 
for some m, n, such that m — n = 0 mod N c . 

Following [2] the normalised chiral condensate can be put in the form 


1 

N s NfN c 


(W) 



(2.5) 


where A (L) is the contributions from all graphs with 2 L links which start and end at some site jc. 
A general graph can be built out of irreducible graphs with less links (if the graph is not already 
irreducible). Specifically, an irreducible graph cannot be separated into smaller graphs which start 

Reducible 

xj 

x 

To obtain the contribution of all general diagrams A(L ) with 2 L links, it is necessary to take 
all possible combinations of irreducible graphs /(/) of 21 links, which form a diagram of 2 L links, 



A(L) = £/(/)A(L-0, L> 1; A(0) = 1, 

i=i 


( 2 . 6 ) 


where the irreducible graphs can begin with an area-0 contribution, a) II, or an area 1 base diagram. 


such as b) —-O , or ... . The first four 1(1) are 


'(!) = 


= I a (l) = 2d, 


I( 2) = 



= / fl (2) = 2J [4(1) ao]. 



+ 



(2.7) 


( 2 . 8 ) 


= 4(3)=2j[4(2)a 0 + 4(l) 2 ^] , 


7(4) = 




+ 2 



+ 
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(2.9) 


= 4(4)+4(4) =2r/[4(3)a 0 + 24(l)4(2)a2+4(l) 3 ^] -4d(d-l)%-, 


( 2 . 10 ) 
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We have defined I a (l) as all iiTeducible graphs of length 21 starting with a) IT, I b (l) as all irreducible 


graphs of length 21 starting with b) 


, etc. The x n are defined as x n = where x n is the number 


of ways of attaching a type x diagram to an area n diagram, defined to reduce over-counting, and 
d x is the total dimensionality of a type x diagram. For example, do = ho = 4^1 l \y More x n 
are defined in appendix A of [7]. In general the 1(1) can thus be put in the form 


1(1) =2dF 0 (l - 1) -4 d(d - l)^Fi(/ - 4) 7 + . 


with 7(0) = 0, 


( 2 . 11 ) 


where F„ (L) represents all possible graphs of length 2 L which start and end on a site on a base 
diagram of area n. The F„ are composed of all possible combinations of iiTeducible graphs which 
add up to 2 L links, 


F n (L)= £ 4,(/ 1 )7 fl (Z 2 )...7 fl (/ / ,)7 fc (k 1 )4(k 2 )...7 fo (^)...flP^..., with F„ (0) = 1. (2.12) 

Zj=l,2,..., 

4=4,8,..., 

Ft ~kj—L 1 


The generating function for all irreducible graphs, including the mass dependence, is 




z=o 
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4m 2 


/(/) =W a + W b +..., 


(2.13) 


where W a is all irreducible graphs starting with an a-type base diagram I, W b is all irreducible 


graphs starting with a fi-type base diagram — , etc. Using (2.12) and (2.11) gives 

2 dx 


W a = 2dx ^ 

n—0 


+ boWb + • • 


W b = -4d(d-l)^j-x 4 


I 

n —0 


ai W a + biW h + ... 


1 -a 0 W a -b 0 W b -.. 

-4d(d-l)^-x 4 


(1-aiWa-bxWb -...) 7 


(2.14) 

(2.15) 

(2.16) 


where jt = — and the “..." contain irreducible graphs starting with higher order (in ] m ) base 
diagrams. The chiral condensate is obtained by taking all possible combinations of all possible 
irreducible diagrams. That is 

JHaL = Um - ( —!— 

N s NfN c 0 m \ 1 — Wj 
It is possible to obtain a simpler system of equations than (2.14) - (2.16) by working in the massless 
limit. One can introduce the variables g x = — such that, taking m —> 0, 


(2.17) 


ga = 


1 

aoga + bogh + ■■■ ’ 


gb = 


Nf 

_)Vc_ 

(a l g a + big b + ...) 7 ’ 


gc = 


Nf 

_ Nc_ _ 

(aiga+b 2 gb + ’ 


.... (2.18) 


The chiral condensate can then be obtained from g = d a g a +d b g b + ..., using 

(W) = 2 

N s N f N c g ■ 


(2.19) 


The prefactors in the numerators of (2.18), and the powers of the quantity in the denominators need 
to be determined for each diagram type. The total contribution of a diagram includes 
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A factor for a number i, of overlapping closed internal loops, 

A mass factor (— 4 ^ 2 )", for n pairs of links, 

(— 1 )* for k permutations of 7 matrices, 

A factor containing the result obtained by performing the group integrations, 
A factor containing the dimensionality of the graph. 


Group integrals for overlapping links of the form 11, or HIT are nonzero V/V r , given by [ 8 , 9, 10, 1 1 ] 

[ dUU a b (Ul) e d = ±-8ffi, (2.20) 

Jsv(n c ) N c 

L m , dvu -- (u,)r '' {u,) » = aw (<'*£+%%) (%• *+) 

+ mhn (<'K - (#,'*£ - %%) ■ 


For finite N c , for example 11 [for SU (3), integrals of the form 


rlTITJ Dl ■■■TJ bNc — ^—P pb v -b Nc 

au u ai u aNc — n ^ ai ... aNc t. 
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ls\J(N c 


( 2 . 22 ) 


are needed. These rules are sufficient to evaluate the diagrams we will use, including 


a ) 


g) 


b ) 


d) 


(2.23) 


□ 


□ 


The specific contributions of these (and other) diagrams are given in [7]. 


3. Group integration with Young Projectors 

To calculate higher order diagrams one needs to evaluate integrals of the general form 


In= / d UUj'...U a /"(U\ S '...(uX 5 " 

JSU(N C ) 


(3.1) 


Any nonzero integral including some combination of U, U 1 can be converted to this form using 


U b } = 


?b\bi % "bNjj a 2 m ' a N 


'ai — \N-\)\ e ‘ 1 iar-a N £ blb2 ‘ ‘ ‘ (U^) b £ and — ^ N \y£aia 2 -a N ^ ^ b 2 " ' “ b N 

Calculating the direct product of n U’s (U f, s) leads to a direct sum of representations R ( S ). The 
integral can be obtained from the Young Projectors P of these representations using [11] 


[ d U R a b (S*) c d = J_(p«)/(pS)/ 8 rs . 

JSU(N C ) d R 


(3.2) 
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Consider for example the integral in (3.1) with n = 2. The direct product N c <g> N c is 



(3.3) 


The Young projectors are thus formed by symmetrising, and antisymmetrising in a,\ and OC 2 , 


pS /3ift _ 1 ( sft xft 

^aioh — u a 2 

The resulting integral is 


I ^(Yi '-'l 


oc\ u a 2 i ? 


pAS pifh }_ l 

^ «l 0.2 2 \ 


sft _ sft c/Jl 

u «l O a2 U ai U a2 


(3.4) 


r _ _pS Si&pS P1P2 _i _p ^5 Si&pAS Pip 2 

2 iV c (iV c + l) “ l0!2 ™ + iV c (iV c -l) ™ ■ 


(3.5) 


More involved examples (I 3 , 14 ) are worked out in [7], 


4. Higher dimensional representations 

Higher dimensional representations can be written in terms of the fundamental and anti¬ 
fundamental. For example, the symmetric (U s ) a b , for a,b = 1, ...,ds, is given by 

(U S )(a ia2 ) [m) = (I ? S ) ai a 2 riYl U Yl Si U Y ^-(F S ) SlS /^ = l - (Ua^Ua^ 1 +U a ^Ua 2 ^) ■ (4-1) 

The antisymmetric (JJ AS ) m n , for m.n = 1,... ,d A s, is given by 

(U A % lCl2l ^- ] = i ( Ua^UaJ 82 - Ua^Ua/') . (4-2) 

The adjoint (U A ) a b , for a,b = 1,... ,dAdj , can be written as 

(t/V = 2Tr(t/f a t/¥) , (4.3) 

where the t u are fundamental generators of SU (N c ) normalised as Tr (t a t b ) = \8 ab . For integrals 
with higher dimensional representation links in the form II, it is sufficient to use 

/ d U (U R ) a b (U m ) c d = -J-8 a d 8 c b . (4.4) 

JSU(N) ClR 


Further considering the adjoint, we are in general interested in integrals with links of the form I I'T , 
for n lines, that is 


I A = / dUU, 


a 1 


‘U a , 


= j dUU ai P'---U a /"Ul S '---Ul S " 


(4.5) 


For example, for n = 3, evaluating the fundamental integral and simplifying using the identity 

CC jTiS ab \si + 1 tlaix'tc T 0 fabdc? results in 


I A - 
y 3 — 


Nr 


-d rfblb-lbi I . 
u a{a. 2 a?, u 1 


_ f 

(Nj - 1)(A? - 4) 1 n c (Af2 - 1 ) /aia2fl3/ 

where if abc = 2 Tr ([t a ,t b ]t c ), d abc = 2 Tr({ t a ,t b }t c ). 


•b\b 2 b-x, 


(4-6) 
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Figure 1: vs. Nf at g = °° including area 1 sub-diagrams up to order J for L = 4,6,8 with 

fermions in the fundamental representation (left), and comparing the fundamental, symmetric, and adjoint, 
including sub-diagrams up to L = 4 (right). 

5. Results 

Results for the normalised chiral condensate are plotted in Figure 1. The solution 

plotted is that which goes to the result of [2, 12] in the Nf —> 0 limit. A more detailed analysis 
of results is presented in [7], A remarkable feature of these results is that as Nf is increased, the 
chiral condensate decreases very slowly and approaches zero as Nf —> °o. Unlike in [3,6], there is 
no indication of discontinuity in any of the solutions obtained. Flowever, we cannot rule out that 
the preferred solution changes at some critical Nf. There are sources of error associated with this 
approach including mis-counting of overlapping diagrams, and over-counting due to symmetries. 
These need to be quantified. For details see [7]. 

References 

[1] J. M. Blairon, R. Brout, F. Englert and J. Greensite, Nucl. Phys. B 180 (1981) 439. 

[2] O. Martin and B. Siu, Phys. Lett. B 131 (1983) 419. 

[3] P. de Forcrand, S. Kim and W. Unger, JHEP 1302 (2013) 051 [arXiv:1208.2148 [hep-lat]]. 

[4] H. Kluberg-Stern, A. Morel and B. Petersson, Nucl. Phys. B 215 (1983) 527. 

[5] P. H. Damgaard, D. Hochberg and N. Kawamoto, Phys. Lett. B 158 (1985) 239. 

[6] E. T. Tomboulis, Phys. Rev. D 87 (2013) 034513 [arXiv:1211.4842 [hep-lat]]. 

[7] A. S. Christensen, J. C. Myers, P. D. Pedersen and J. Rosseel, arXiv: 1410.0541 [hep-lat]. 

[8] K. G. Wilson, “Quarks and Strings on a Lattice,” CLNS-321. 

[9] I. Bars and F. Green, Phys. Rev. D 20 (1979) 3311. 

[10] M. Creutz, Cambridge, Uk: Univ. Pr. ( 1983) 169 P. ( Cambridge Monographs On Mathematical 
Physics) 

[11] P. Cvitanovic, Princeton, USA: Univ. Pr. (2008) 273 p 

[12] P. de Forcrand and S. Kim, Phys. Lett. B 645 (2007) 339 [hep-lat/0608012]. 


7 







